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We study the phenomenology of a Z'-boson field coupled to hypercharge. The Z' propagator has 
a nontrivial Kallen-Lehmann spectral density due to the mixing with a higher-dimensional inert 
vector field. As a consequence detection possibilities at hadron colliders are reduced. We determine 
the range of parameters where this field can be studied at the Tevatron and the LHC through its 
production cross section via the Drell-Yan mechanism. 
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O ■ I. INTRODUCTION 

o ■ 

^ I The present day data from high-energy colhders hke LEP and Tevatron show that almost all measurements are 
Q , described by the standard model at the loop level. Therefore extensions of the standard model tend to be strongly 
constrained. Typically such extensions will spoil the agreement with the data through a variety of effects, one of the 
most important of which is the appearance of flavor changing neutral currents. Even the most popular extension, the 
^ minimal supersymmetric extension of the standard model, has to finely tune a number of parameters. This leaves only 

one type of extensions that are safe, namely, extensions with singlet particles, since singlets do not affect precision 
■ variables through loop effects. Since the discovery that neutrinos are massive, it is clear that singlet fermions play a 
Qh| role in nature. Moreover arguments from cosmology suggest that singlet scalars could be important [H, H, S, 0, 11| • In 
f-^ • a recent analysis of the Higgs-search data from LEP2 it has been pointed out that such singlet scalars may already 
^ ' have been seen as a smeared-out Higgs boson [6] . Given this situation it is natural to ask whether also singlet vector 
, bosons can be present. 

, Therefore, we decide to study renormalizable extensions of the standard model containing extra gauge bosons, but 
^ ' no extra fermions or scalars. Demanding the absence of anomalies in the gauge currents, additional vector bosons can 
\ couple to linear combinations of hypercharge (Y) and the difference of baryon and lepton quantum numbers (B-L). 
OO i However, they can have an arbitrary mixing with the standard model hypercharge field. 

In this paper we limit ourselves to the case where the vector bosons couple to the hypercharge only. Such models 
" [ , are, after a redefinition of the gauge fields, equivalent to the standard model with a number of inert gauge bosons, 
00 ' that only couple to the standard model via the mixing with the hypercharge field. Since the extra vector bosons 
\ are singlets, their mass does not necessarily have to arise from spontaneous symmetry breaking, but can be put in 
by hand without violating renormalizabilitv. For phenomenological purposes the differences are minor. Both cases 
have been studied in the literature 0, [li H Elj [13, [13, [3 • However if the mass is not generated by the Higgs 
mechanism, an interesting possibility arises when an infinity of fields is present. It is possible that the extra vector 
bosons move in more than four dimensions [l5| . This will not violate the principle of renormalizability, if the number 
of dimensions is smaller than or equal to six; also fractional dimensions are possible. As a consequence one finds, 
d ' after diagonalizing the mass matrix, a physical neutral Z'-boson field coupled to hypercharge. However, this field 
does not have a single particle peak but has a well-defined nontrivial continuous Kallen-Lehmann spectral density. 
Therefore, the search for such a field at hadron colliders becomes more complicated, since the signal is smeared out 
over a large invariant-mass range and no clear peak is present. Because of their higher-dimensional nature we call 
this class of models Heidi models (high-D for higher dimension). 

It is the purpose of this paper to determine the range of parameters of the model that can be studied at hadron 
colliders. We consider a minimal version of Heidi models with a reduced set of free parameters, namely, an extension 
of the standard model with one single additional inert gauge boson mixing with the standard model hypercharge 
field, which affects the usual standard model neutral gauge bosons. While the photon is protected from further 
mixing by the exact electromagnetic symmetry, no conserved quantum number prevents the Z boson to mix with the 
additional Z' boson. Theoretical calculations of various electroweak observables are then modified by the presence 
of the additional field, which may change the fits of the electroweak precision data and affect the bounds on the 
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Higgs mass [T^ [13, [IB [13 ■ According to this point of view, measurements from LEPl and SLD at the Z pole 
[l7j were analyzed in Ref. [l^l, both for a four-dimensional and for the Heidi Z' models. The analysis on the 
technical level is not significantly different for the four-dimensional and the higher-dimensional case. The results of 



this analysis were recently confirmed in Ref. Fitting the data leads to constraints on the allowed parameter 



range of the Heidi models. We will restrict ourselves to the limits given in Ref. p^], which are on the conservative side. 

The Drell-Yan-like production of a lepton pair at hadron colliders plays an important role in the present and future 
experimental program at the Fermilab Tevatron and CERN LHC colliders thanks to its large cross section and clean 
signature of the final state. This is also the primary discovery mode for additional neutral gauge bosons and has 
already been widely studied in the literature for the Tevatron and the LHC [3 [l^, 113, HHj uf\ . It has also been 
shown that the forward-backward asymmetry due to the interference between the three different exchange modes (7, 
Z, and Z') could be an interesting observable [1^ HI, H^. Therefore, in this work, we investigate for the first 
time the possible effects of a Heidi Z' on the Drell-Yan dilepton invariant-mass spectrum and the forward-backward 
asymmetry, both at the Tevatron and the LHC. 

Because of the important theoretical uncertainties coming from QCD effects and missing large higher-order 
contributions, leading-order results for the Drell-Yan process, with or without additional Z', are quite unreliable. 
A systematic approach to this problem is based on perturbation theory truncated at the next-to-leading order 
(NLO) [13, mi or next-to-next-to- leading order in the strong coupling constant a^. In kincmatical regions where 
higher-order contributions are enhanced due to the soft and coUinear parton emission, soft-gluon resummation to all 
orders in can be performed and matched with the fixed order predictions in order to improve the description of 
observables under consideration. However, a general procedure has not yet been developed and the resummation has 
to be done explicitly process by process. An alternative way of treating the soft-gluon emission is the parton shower, 
which approximates the full resummation calculation. However, the parton shower is universal for all processes 
and can be implemented easily in the analysis chains of the Tevatron and LHC experiments. Consequently, we 
implement in this work the Heidi Z' bosons in the Monte Carlo generator MC@NLO [30], generalizing the modified 
MC@NLO program including Z' gauge bosons inspired by grand unified theories of Ref. [20*1, which allows us to 
match NLO perturbative calculation with the parton shower of the Monte Carlo generator HERWIG 31] including 
color coherence effects and azimutal correlations within and between jets. We subsequently compare the Monte Carlo 
predictions to the experimental results from the Tevatron for the invariant-mass spectrum and forward-backward 
asymmetry of the Drell-Yan production of lepton-pairs. In addition we study the sensitivity of the LHC experiments 
to the Heidi models for this channel. 

The remaining part of the paper is organized as follows. In Sec.|TT]we first describe the theoretical framework, i.e. the 
minimal Heidi Z' model, and discuss the implementation of Z'-boson production in the MC@NLO generator. Section 
mil is devoted to the numerical analysis of the experimentally allowed Heidi parameter space and the investigation of 
the sensitivity of the LHC experiments to the model. Our conclusions are given in Sec. IIVI 

II. THEORETICAL FRAMEWORK 
A. The minimal Heidi model 



There is a large class of models containing extra neutral vector bosons [21|,|24]. The simplest extension consists of the 
standard model plus one additional C/'(l) symmetry, the extended gauge group being S'C/(3)cxS't/(2)iX[/(l)y xJ7'(l). 

We associate with the extra symmetry a Z'^ field that is a singlet under the standard model gauge group. Moreover, 
we assume that no particle carries a V (X) charge. Therefore, this Z'^ field is inert and its presence can only be noticed, 
because it can mix with the Abelian field -B^ associated with the U{X)y symmetry. Because it does not couple to any 
fermions or Higgs fields, the additional Z'^^ field is allowed to move in D dimensions. As long as D < 6 the theory 
will stay renormalizable. In contrast, the field must be four-dimensional, because of its couplings to the fermionic 
fields. Therefore, in this basis, the hypercharge basis, we have two fields, a four-dimensional field coupled to the 
hypercharge and a higher-dimensional field coupled to nothing. However, the two fields can be mixed by a mass term. 
In order to describe the physics correctly, one has to transform to the mass basis. Allowing the ?7(1) fields to have 
masses and mass- mixing terms, one finds for the gauge field part of the Lagrangian in the hypercharge basis [l^ : 

- ll„ = Ib^^B^-^ + \mlB,B>^ + \z[^,Z'^^ + \mlZ'^Z'^^ + M^^^ Z'^B^^ . (1) 
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We have a four-dimensional mass TO4, a higher-dimensional mass ttid, and a scale Mmix that connects the four- 
dimensional fields with the higher-dimensional ones. The two massive hypercharge fields Z'^ and are transformed 
to the mass-eigenstates Z'^^ and _B^, the latter being massless since the electromagnetic symmetry is exact. As a 
consequence, the three parameters 7714, mc, and Mmix are not independent and must satisfy the condition 



for Z) 7^ 6 
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where we have introduced the quantity /iihd in reference to the mixing of the low-dimensional and high-dimensional 
fields 
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By compactifying the higher dimensions and subsequently taking the continuum limit, one can derive the 
hypercharge-boson propagator which is of the form 
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and its nontrivial corresponding Kallen-Lehmann spectral density 
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which becomes in the four-, five- and six-dimensional cases 
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respectively. For the three choices of the number of dimensions D, we recover models with one massless excitation 
which becomes the photon after the breaking of the electroweak symmetry. This is guaranteed by the condition of 
Eq. ([2). However, while in the four-dimensional case we have a second peak corresponding to a massive resonance 
located at 
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for five- and six-dimensional singlet fields, we now have a massive continuum. As stated above, the presence of this 
new field affects the standard model Z-boson through a mixing which can be described by a quantity ay defined by 
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for the four-, five-, and six-dimensional cases, respectively. This quantity enters directly the vector and axial-vector 
coupling strengths of the Z boson to fcrmions, which read now 

Vj = e/(sin^ - ay) - - ay) , 

«/ = ■ (9) 

and e/ denote the weak isospin quantum number and the electric charge of the fermion /, and svy? Qw is the 
squared sine of the electroweak mixing angle. The coupling strengths of the new Z' state to fermion are simply given 
by the hypercharge 

a!j = F(m4, m£)) sin^H'-^, (10) 

with a prefactor F depending on the four-dimensional and high-dimensional masses. 

One can interpret the propagator of Eq.(4) in two ways. The simplest way to treat the theory is to consider the field 
as a hypercharge gauge field with a nontrivial Kallen-Lehmann spectral density. The other interpretation is to say 
that one has many (in this case a continuum of) single hypercharge coupled Z' fields, however, with a reduced coupling 
constant. The latter interpretation would be natural if we had an integer number of compact higher dimensions. 



B. Heidi ^'production at hadron colliders 



The Drell-Yan-like production of a lepton pair plays an important role in the present and future experimental 
program of hadron colliders thanks to its large cross section and clean signature of the final state. Since any observables 
related to this process are highly sensitive to the existence of an additional neutral gauge boson, this is their primary 
discovery mode. Following the conventions for the Drell-Yan process of Rcf. [2g], which are the ones used in the Monte 
Carlo generator MC@NLO, the corresponding squared matrix element reads 



qq or qg 
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(11) 



and can be written as 
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where M is the invariant-mass of the produced lepton pair, Ci is the corresponding color factor and the functions 

\A,,Bi 

the coefficients 



Ti \ a'' a depend on the kinematics of the process and on the order of the perturbative calculation. We have defined 



4') = a^p'+vf\ 



(13) 
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FIG. 1: The (m4,mo) plane for Heidi models with a four-dimensional (left), five-dimensional (center), and six-dimensional 
(right) Z' boson. We show the regions excluded by the LEP collider with (light gray region) and without (below the black line) 
leaving the forward-backward asymmetry of the bottom quarks, by the Tevatron (black region). We present also the region 
(dark gray) that the LHC experiments will be able to exclude at the 95% confidence level with an integrated luminosity of 100 
fb-\ 



and a modified propagator obtained by removing the photon and Z contributions from tlie propagator of 

Eq. (dl, since they have already been taken into account. The squared Z'-boson exchange as well as its interferences 
with the photon and Z-boson exchanges are included. We have adapted the MC@NLO program of Ref. [20j| for the 
Drell-Yan process including Z' inspired by grand unified theories to the Heidi model case, which allow us to match 
the matrix elements of E g. (IT^ . describing correctly the hard parton emission by the initial state, to the HERWIG 
parton shower algorithm |3l| describing the soft and coUinear emission. 

In the considered process, the incoming quark, antiquark, and gluon can give rise to an initial-state parton shower 
which is modeled in HERWIG by starting with the hard scattering partons and reconstructing backwards the preceding 
branchings. For a specific splitting of partons i — » jTc, the energy fraction of parton j is distributed according to the 
LO DGLAP splitting functions and the phase space is restricted according to an angular ordered emission which is 
based on the Sudakov form factor. This allows us to sum the virtual corrections and unresolved real emissions to 
all orders and to correctly treat the infrared singularities. The parton shower stops when a cutoff scale defined by 
(2.5GeV)^/i?f is reached, where a nonperturbative stage is imposed. 

In order not to double count any contribution, the matching of the NLO matrix element to the parton shower 
achieved in MC@NLO uses two separate samples of Born-like and hard emission events which can have positive or 
negative weight. They are statistically distributed according to the positive-definite corresponding part of the NLO 
cross section and made separately finite by adding and subtracting the NLO part of the expanded Sudakov form 
factor fs^l. The total cross section is then given by their weighted sum averaged over the total number of events. 



III. ANALYSIS 



A. Experimental constraints 

Although there is no experimental evidence for the existence of a Z' boson, a number of data can be used to 
constrain the parameter space. For example, limits on the parameter ay, and thus indirectly on 7714 and tod, can be 
obtained from precision measurements at the Z pole, strongly constraining the mixing between the standard model 
Z boson and the new field. 

Regarding the limits from the LEP experiment there are two ways to proceed. One can take all the data from 
the precision measurements |17i] and try to make a fit. The problem here is that the standard model is only barely 
compatible with these data due to the inconsistency of the forward-backward asymmetry of the bottom quarks 
with the rest of the data. Adding a Z' boson does not improve the situation, one finds a limit ay < 1.8 10~* at the 
95% confidence level. This problem with the data has led a number of authors to consider fitting the data leaving 
out Apg [H, [3, [H, m, m, m, H^l. in this case a good fit to the data is possible if a Z' boson is present p^. [l3|. 
With a Higgs mass ranging from 115 to 495 GeV, a larger range of ay is allowed: ay < 12 10"''. As the situation is 
somewhat controversial we will present both limits in the following figures and discussions. 
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TABLE L Heidi benchmark points accessible at the Tevatron and aUowed by the LEP constraints after leaving out the forward- 
backward asymmetry of the bottom quarks. 
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Imposing these limits on minimal Heidi models with three free parameters, the number of dimensions D, the 
four-dimensional mass and the high-dimensional mass rn_D, leads to the excluded regions shown in Fig. [T]for a 
four-dimensional (left), five-dimensional (center), and six-dimensional (right) additional neutral gauge boson. The 
light gray regions are the regions where ay > 12 10^^, while the dark line comes from the more stringent limits 
obtained by keeping all LEP data, i.e. it corresponds to ay — 1.810-4. The excluded regions are thus those below 
that line. 

The best limits on the presence of a four-dimensional Z' boson come from the Tevatron collider. The principle 
guiding the search is straightforward. One uses the decay of the Z' boson into an electron-positron pair and looks 
for a peak in the invariant-mass. In addition one can use information on the forward-backward asymmetry of the 
leptons. One therefore needs a prediction for the Drell-Yan production cross section. Moreover, when one is working 
within a specific model, one can use the distribution of the leptons in the center-of-mass angle cos 0* with respect to 
the beam axis, leading to somewhat different limits for different models. This search has been made at the Tevatron 
for the class of models discussed in ReL 18], but not for the specific models presented here. 

The model dependence is contained in the values of the vector and axial- vector coupling constants. As argued in 
Ref. 18] it would be useful to have lower limits on the mass of the Z' boson presented in the coupling constant plane. 
Unfortunately such a comparison has not been presented in the literature. Instead we estimate the lower limits for 
the allowed mass of the Z'-boson as a function 7714 and nijj , as follows. We use our implementation of the Heidi model 
in MC@NLO to predict the total cross section as a function of TO4 and niD- Then, this cross section is normalized 
to the sequential Z' boson, which has a lower mass limit of mz' = 923 GeV [35|. We then connect the lower bounds 
on mz' within the two studied models. The results are given in the dark area of Fig. [T] (left panel). The derived 
bounds are of course somewhat qualitative, because the sensitivity to the cos 0* distribution is not exactly modeled 
this way. A precise analysis would require taking into account a bidimensional distribution, including the angular one 
as well, whereby one cannot ignore detector effects [11]. However, this needs a detailed simulation of the detector and 
comparison with the actual data, which is beyond the scope of this paper. 

For the higher-dimensional cases this procedure is not sufhcient and a bin-by-bin comparison of predicted and 
measured Drell-Yan cross sections is needed. Unfortunately, these measurements are only presented in the literature 
[STj for a limited integrated luminosity. Within this limited statistics the Tevatron gives no improvement over the 
LEP data. However, we note that in Ref. 36] there appears to be an excess in the 300 — 400 GeV range, which might 
be consistent with a spread-out Z' boson. We hope that published numbers for the Drell-Yan cross section with a 
larger integrated luminosity will appear soon, so that these questions can be studied in more detail. 

In Fig. [21 we compare our results with the limited statistics data published in Ref. [13] ■ We show both the differential 
cross section da/dM and the forward-backward asymmetry. As an example, we chose the two Heidi benchmark points 
presented in Table IH which are just within the region of the parameter space allowed by the constraints coming from 
LEP, provided that the forward-backward asymmetry of the bottom quarks is left out, and "friendly" for the Tevatron. 
We show our predictions for the Heidi model together with those of the standard model, i.e. when no additional Z' 
boson is present. We use the CTEQ6M (NLO MS) ^] sets for the parton densities in the proton and antiproton. 
As can be seen from the figures, all the predictions are compatible with the data for both observables. However, data 
with a higher luminosity would make it sensible to divide the studied invariant-mass range in smaller bins, which 
could be used to emphasize deviations from the standard model. 

B. Predictions for the LHC 

In this section, we discuss the sensitivity of the LHC experiments to Heidi models, considering four-dimensional, 
five-, and six-dimensional extra neutral gauge bosons. The discovery potential of the LHC is determined through the 
investigation of possible deviations from the standard model predictions in the tail of the Drell-Yan invariant-mass 
distribution, the studied range being 
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FIG. 2: Drell-Yan differential cross section and forward-backward asymmetry with and without Heidi Z' contributions at the 
Tevatron. The results are presented for the two benchmark points of Table H] 



TABLE II: Typical Heidi benchmark points accessible at the LHC and allowed both by the Tevatron and by the LEP constraints, 
after leaving out the forward-backward asymmetry of the bottom quarks. 
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The main source of background for Drell-Yan l^l^ events is QCD multijet and direct photon production where 
the jets have a large electromagnetic component. However, with typical experimental selection cuts one can obtain a 
signal almost free from any QCD background, its part to the measurements being reduced to at most 1% [39l[40l[4]|. 
For that reason, we apply the following cuts 

j PT > 20 GeV 
[ \v\ < 2.5 

to both leptons, which allows us to consider as a new physics signal any significant excess of events with respect 
to the standard model DrcU-Yan irreducible background, neglecting thus any other source of background. 

For each bin i of the investigated range of Eq. p^ . we calculate the quantities and Ni, representing the 
number of events related to Drell-Yan-like lepton-pair production with and without an additional neutral gauge 
boson, respectively. For a given luminosity L, they can be derived from the differential cross section by 

where Mi (M2) is the lower (upper) limit of the considered bin, and da {da') the Drell-Yan cross section without 
(with) an additional Z' gauge boson. Expecting the standard model Drell-Yan background, we can estimate the 
significance of a possible Heidi signal through the quantity 
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FIG. 3: Drell-Yan difFerential cross sections with (solid red) and without (dashed blue) additional Heidi Z' bosons at the LHC 
for the benchmark points presented in Table HIl The typical experimental cuts of Eq. (|15p are applied. 



where n is the number of bins which the range of Eq. p4|) is divided in. The LHC will then be able to exclude a 
given point of the Heidi parameter space at the 95% confidence level if the corresponding is bigger than 21. 

In order to take advantage of the differences between the standard model and the Heidi predicted line shape of the 
Drell-Yan mass spectrum in the high invariant-mass regions, we divide the studied range into nine bins of 100 GeV 
for 200 GeV < M < 1100 GeV and we take three larger bins at very high invariant-mass in order to have higher 
statistics, i.e. 1100 GeV < M < 1300 GeV, 1300 GeV < M < 1600 GeV, and 1600 GeV < M < 4500 GeV. That 
allows us to avoid spurious sensitivities since we have then at least 50 expected events for each bin. This way we 
have a good idea of the typical reach of the LHC in the parameter space of the models. The results shown in Fig. [1] 
use cross sections calculated at a center-of-mass energy of 14 TeV with again the CTEQ6M (NLO MS) [111 sets of 
parton densities, assuming a LHC lumisotiy of 100 fb~^. 

For the four-dimensional model, the LHC will be able to reach a mass range of several TeV, as can be seen from 
the dark gray region in the left panel of Fig. [T] This range can be divided in two areas, one at small value of 
TO4 and one at larger values. The first zone corresponds to a region where the Z' mass is smaller than or equal 
to about 3 TeV and is fairly typical since it is also found for other models with electroweak size couplings. In 
contrast, the LHC-reachable zone where m4 is larger corresponds to a heavier Z' but with a much larger width 
due to the large value of the mixing parameter ay, which makes it possible to detect the new vector boson 
via its interferences with the standard model photon and Z boson. In comparison to the results of previous 
experiments, the LHC will allow us to considerably enlarge the model parameter space that can be investigated in 
the four-dimensional case. For the five-, and six-dimensional cases, however, the LHC will not be able to cover a 
parameter space region much larger than is already constrained by the LEP collider. This is the case if we take 
the strict limits that exist if we keep the data coming from the forward-backward asymmetry of the bottom quarks 
in the electroweak precision fit. As stated before, the resonance is replaced by a massive continuum, rendering the 
detection of any excess in the Drell-Yan distributions more complicated and restricting the regions of the parameter 
space that can be reached at the LHC, as can be seen in the central and right panels of Fig. [TJ However, if one 
prefers to take the more conservative limits obtained by leaving out the forward-backward asymmetry of the bot- 
tom quarks, the LHC is sensitive to a range of parameters roughly double the size that is reached by previous colliders. 

As an example, we show in Fig. [3] lepton-pair invariant-mass spectra for the benchmark points presented in Table 
im We compare MC@NLO predictions including a Heidi Z' boson (solid red curves) with the ones expected in the 
standard model case (dashed blue curves). In the four-dimensional model (left panel of Fig. [3]), far from the peak 
region located at the Z' mass given by Eq. ([7]), the spectra with and without Z' will of course coincide, while around 
the Z' mass, we have a clear resonance due to the presence of the extra field. For five-dimensional (central panel of 
Fig. [3]) and six-dimensional (right panel of Fig. [3]) extra neutral gauge bosons, the resonance is replaced by a massive 
continuum, which makes the detection a bit more tricky. However, large deviations can still easily be detected if we 
consider the line shapes of the distributions rather than the total number of events. 
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IV. CONCLUSIONS 



We studied a new class of renormalizablc models containing a Z' boson with a nontrivial Kallen-Lehmann repre- 
sentation, derived from the mixing with an inert higher-dimensional field. Existing LEP limits constrain this class 
of models, but the amount depends on the interpretation of the analysis of the precision data, allowing for strict or 
loose limits. In the case of the loose limits we found, using a state of the art MC@NLO program, that the LHC will 
be able to study a large range of the allowed parameter space. In the case of the strict limits even the LHC will have 
difficulty to improve on existing bounds, but there may be some hope for a large luminosity option. Present bounds 
from the Tevatron, which may be relevant at low masses, are not satisfactory, since only a small fraction of the data 
have been analyzed for the Drell-Yan cross section. We hope the situation will improve soon. 
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